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1. INTRODUCTION 


Electrical power is one of the most important infrastructure components for the economic growth 
and welfare of the developed countries. The demand for electrical power has increased rapidly due to several 
reasons [1]. In order to meet the increasing power demand, massive addition to the installed generating capacity 
is required. Due to the unique characteristics of the electricity, independent system operator (ISO) and regional 
transmission organization (RTO) are responsible to balance the power generation and demand [2]. In real- 
time power system operation, ISO forecasts and schedules generation to assure that the sufficient generation 
and back-up power is available to meet unexpected demand or generation loss. It must be a non-commercial 
organization, neutral and independent of commercial players. Currently, Nine ISO/RTO are responsible for the 
operation in the electricity markets. 


In the deregulated power markets, the bids are submitted by the Genco’s to the power exchange (PX) 
to buy and sell of the electric power. They gradually build their offers strategically to intensify their profits. 
This process is called bidding strategy [3]. The electricity market has been ruled by the bidding strategies. 
These strategies also grab the bidder’s attention with the advantage of increasing profits. Flow of information 
in the electricity market is given in Figure 1 [4]. 
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Figure 1. Flow of information in the electricity market 


Several researchers in the past five decades used various mathematical techniques such as conven- 
tional, naturally inspired algorithms and hybrid methods for solving bidding strategy in the electricity markets. 
Highlights of these techniques are provided in [5]. The emerging electricity market has been analyzed by ap- 
plying game theory in [6]. Network optimization technique is adopted for Nash equilibrium bidding strategy in 
[7]. The bidding strategy problem has been solved effectively using Monte Carlo simulation [8]. Discrete-state 
and discrete-time Markov decision process has been applied for optimal multi-period bidding strategy [9]. A 
Lagrangian relaxation [10] based approach has been adopted for strategic bidding. Strategic bidding problem 
for competitive power suppliers in the England - wale electricity markets is provided in [11, 12]. These al- 
gorithms such as genetic algorithm (GA) [13], particle swarm optimization algorithm (PSO) [14], differential 
evolution (DE) [15], invasive weed optimization (IWO) [16] and krill herd algorithm (KHA) [17], bacterial 
foraging algorithm [18], agent-based algorithm [19], bat inspired algorithm [20] have been applied in solving 
the DSOBS problem. Advantages and limitations of these methods are explained by few researchers. Some 
other algorithms [21, 22] by combining two or more algorithms have been developed to get global optimal 
solution for the DSOBS problem due to the complexity involved in the problem. 

Bacterial foraging optimization (BFO) algorithm has been applied to solve various complex engineer- 
ing problems. Itis inspired by the social foraging behavior of Escherichia coli. Biology behind the foraging 
strategy of E. coli is used as an optimization algorithm. Scholars have been using the BFO algorithm with other 
available methods to enhance the local and global search properties. The BFO algorithm has been applied in 
solving several real time problems in engineering. Also, it has been found that the BFO algorithm is capable to 
provide the global solution. In contrast, mathematical modelling and modification of the algorithm is a major 
part of the research on BFO algorithm. Better solution can be achieved by the proper selection of the control 
parameters and updation of the algorithm. These observations motivated to introduce the BFO algorithm for 
solving single sided optimal bidding strategy in this article. 

The contributions of the article are listed below : 1) Bacterial foraging optimization algorithm has 
been proposed for solving suplier side optimal bidding strategy. 2) The proposed algorithm has been coded in 
MATLAB (2019 version). 3) The BFO algorithm is tested on IEFE 30 bus system and results are provided. 4) 
The results like powers, profits and marginal cost price of the proposed method have been compared with the 
existing techniques. 

The proposed bacterial foraging optimization algorithm is developed in MATLAB (Version 2019) and 
tested on standard test case available in the literature. The code has been executed on personal laptop (8 gb 
RAM, intel i5 processor, 2.3 GHz). The remaining paper is organized as follows: section 2 is bidding strategy 
problem formulation, section 3 is bacterial foraging optimization, section 4 is development of the algorithm, 
section 5 is case study, and section 6 is conclusion. 


2. PROBLEM FORMULATION 
In this section, a brief description about supplier side bidding strategy is provided. Also, mathematical 
formulation is given. The above are described in the following sections. 


Supplier side optimal bidding strategy for... (Uday Kiran L) 


60 o ISSN: 2252-8792 


2.1. Bidding strategy 


In an open electricity market, the hourly aggregate supply is matched to hourly market clearing price 
(MCP) by the market operator. An optimal bidding strategy is mandatory for the Genco’s in order to maximize 


its profits. In the deregulated power markets, Genco’s build their offers to maximize profits. It is called strategic 
bidding. 


2.2. Mathematical formulation 

At the equilibrium point, The ISO matches the price for suppliers and consumers. At this point, the 
price is same for both the players. This price is called marginal cost price. Bidding curve of a suppliers is 
provided in Figure 2. The curve is mathematically modelled as follows: 


ai + b;P; = MCP i=1 ton (1) 
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Figure 2. Bidding curve of a suppliers 


2.2.1. Objective function 
The profit of supplier and consumer is determined using the following equations. 


Supplier; profit = MCP x P; — C;(P;) (2) 


here, supplier; cost function is 
Ci(Pi) = eP; + fiP? (3) 


In an ideal case, the power generation is equal to the power demand. In the pool market, when the 
price elasticity is present, the power demand varies. The mathematical equation of the cumulative fore-casted 
pool demand (QMCP) is 


Qumor = Qo — K x MCP (4) 
2.2.2. Equality constraint 
The power equality constraint is 
XOP. = Quer (5) 
i=1 
Expression of the marginal cost price (MCP) is derived from the above equations and provided below. 
Qo + Sat b 
IO an 


The expressions for power and load are derived in terms of the MCP and bidding co-efficients and 
given blow. 


MCP = (6) 


p, = Yra (7) 
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2.2.3. Inequality constraints 
Generator inequality constraint is 


2.3. Probability density function (pdf) 

Consumers and suppliers set the MCP to maximize their profits when they know their bidding coef- 
ficients. Each GENCO can guess their competitors bidding coefficients using the probability density function 
(pdf). The expression for the pdf is given below. 


1 —1 
dfp(Ti, Yi) = ! 
Meo) ran A E A k 
here, 
OBUP a uP a alos — ot Ilo = a) ig 
aur = Fey) toy | nA) oR 


It is observed that the bidding co-efficients cannot be selected directly. (i) Here, first objective is to 
maximize the pdf by selecting the bidding co-efficients. (ii) The second objective is to maximize the profit 
of the suppliers and consumers. The above two observations are motivated to introduce two stage bacterial 
foraging algorithm. The description of the BFO algorithm and development of the proposed approach are 
given in the proceeding sections. 


3. BACTERIAL FORAGING OPTIMIZATION (BFO) ALGORITHM 

BFO [23, 24] algorithm is one of the biologically inspired optimization algorithms. In this section, 
the details such as background of bacteria and its foraging technique, short description of development of BFO 
algorithm and implementation of the said algorithm are provided. 


3.1. Background of bacteria and its foraging technique 

Bacteria chemical factories capable of bringing about significant changes in nature. BFO algorithm 
is based the foraging techniques of the E. coli bacterium cells present in the intestine of the human digestive 
system [25, 26]. In the race of foraging for food location and survival, these cells with poor foraging techniques 
will be eliminated and those with good foraging techniques obtain enough food and reproduce their genes. As 
this process gradually intensifies the bacterium cells with poor foraging techniques are either extinct or re- 
designed. These foraging techniques and evolutionary principles helped scientists to hypothesize and correlate 
them with the optimization processes. A foraging animal always strives to ingest and gain energy by overcom- 
ing the physiological and environmental defects like poor sensing, cognitive capabilities and high density of 
prey, risky predators respectively. 


3.2. Development of BFO algorithm 
Typically, the BFOA consists of four main mechanisms: chemotaxis, swarming, reproduction and 
elimination-dispersal event. 


3.2.1. Chemotaxis 

In this process the movement of bacteria is enabled by helix-shaped flagellum. These set of rigid 
flagella helps the bacteria to swim. These flagella help the bacterium cells to move either by pushing them 
while rotating in an anti-clockwise direction by pulling them while rotating in the clockwise direction. This 
mechanism which enables the flagellum to spin by creating rotational forces is called bacterium motor. The 
biological motor can make the bacterium cells to swim or tumble. The E. coli bacterium cells maintain a spec- 
ified direction during swimming and a random direction during tumbling. These alternate modes of operation 
are carried throughout its lifetime which enables them to search for nutrients. 


Ati) 


(j +1,k, 1) = 6 (j, k,l) + C(i) x TAT@AD 


(1) 
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3.2.2. Swarming 


In this process, the bacterium cells swarm together just laid a group of swarming birds flying towards 
a certain direction during the migration season. When the bacterium cells reach a suitable location they start 
attacking other cells and converge at that particular location. This swarming pattern forms based on the dom- 
inance of the two stimuli relative distance between the respective bacterium from the signalling bacterium the 
performance of swarming from the signalling bacterium. The performance of swarming is to create a unique 
high dense bacterium for better foraging. The concentric patterns are a result of merged bacterium cells at the 
fittest point. 


Jecl (j, k, 1), PG, k,1)) => BAO RL), PG, k,1)) (12) 
: 2 P 
Jecl’ (j, k, l), PU, k,1)) = So [-dattract x egp attract] X 0” T al (13) 
i=1 m=1 
s p 
st S“[-Arepeitent x erp “renettent] X [0™ z or]? (14) 
i=1 m=1 


3.2.3. Reproduction 

It is applied after chemotaxis, the process is provided below: (i) Fitness values of bacteria are store in 
ascending order, (ii) Bacteria with least fitness value are eliminated, (iii) The remaining are divided into two 
identical and filled at the same location. 


3.2.4. Elimination and dispersal 

The Process of swarming is affected when a group of bacteria gradually or suddenly reaches a new 
location due to the consumption of nutrients or by any other influence like attracting to fittest bacterium. This 
results in elimination of set of bacteria or dispersion from the respective group. The flowchart of the BFOA is 
provided in Figure 3. 


Figure 3. Flowchart of the BFOA 
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4. DEVELOPMENT OF THE BFO ALGORITHM 

The problem of supplier side optimal bidding strategy (DSOBS) is formulated in section 3 in this ar- 
ticle. The objective is maximization of the profit of suppliers by selecting the bidding coefficients strategically. 
It is found from the survey of literature that the profit cannot be maximized by choosing all bidding coefficients 
independently. One of the bidding coefficients independently can be fixed. Now, the other bidding coefficient 
is determined by using any optimization method. This article presents the development of BFO algorithm for 
maximizing the profit of suppliers and consumers by choosing the appropriate bidding coefficients. 
Step 1: ais kept constant and b is initialized. Range of the value of ’b’ is [b 10b] 
Step 2: The objective function is maximization of ‘profit. Again, BFO algorithm is used with control parame- 
ters shown in Table 1. 
Step 3: Calculate 
[1] MCP using (6) 
[2] Power using (7) 
[3] Profit of supplier using (2) 
Step 4: Check equality constraint using eqn. (5), if ‘yes’ go to ‘endelse. 
Step 5: Apply operators of BFOA 
Chemotaxis [Refer equation 12] 
Swim [Refer equation 17] 
Reproduction [Refer section 3.2.3] 
Elimination [Refer section 3.2.4] 
Step 6: Update the values of ‘b’ go to step 3. Repeat the steps from step 3 to 5 until either the maximum 
iterations are reached or equality constraint is met. 
Step 7: Final result. 


Table 1. Control parameters 


SNo Control parameter Value 
1 Population size (S) 20 

2 No. of Chemotaxis steps (Nc) 30 

3 Ns 4 

4 No. of repro steps (Nre) 10 

5 No. of elim-disp (Ned) 5 

6 cá) 0.05 
7 Probability of elimi-disp (Ped) 0.02 


5. CASE STUDY 

The of the proposed algorithm has been developed in MATLAB (Version 2019). Applicability and 
simulation results of the proposed algorithm in comparison with existing algorithms has been done. Effective- 
ness of the proposed BFOA is tested on a standard system. Data of the system is available in [13] and presented 
in Table 2. The code has been executed on Personal laptop (8 GB RAM). The control parameters and their 
numerical values are shown in Table 1. 


Table 2. Generator data 


Supplier e f Piin Porai 
1 2.00 0.00375 20 160 
2 1.75 0.0175 15 150 
3 1.00 0.0625 10 120 
4 3.25 0.00834 10 100 
5 3.00 0.025 10 130 
6 3.00 0.025 10 130 


5.1. Case 01 

BFO algorithm is tested on the IEEE 30 bus system. In this case, the values of Qo is 500 MW and K 
is 0. During the execution of the algorithm, the code is run for 20 times. Bidding coefficients, output powers 
and profits of the case study are presented in Table 3. 


Supplier side optimal bidding strategy for... (Uday Kiran L) 


64 o ISSN: 2252-8792 


Table 3. Simulations results of IEFE 30 bus system for case 01 


Supplier a b Powers Profit 
1 2.00 0.038452 160 1208.351 
2 1.75 0.144258 58.24431 430.013 
3 1.00 0.211042 43.36678 279.3588 
4 3.25 0.082114 84.05669 521.2492 
5 3.00 0.12672 56.44106 324.0376 
6 3.00 0.155882 45.88207 275.5284 
MCP 10.15219 
Profit 3038.538 


It is observed from Table 3 that 
— All output powers are within generator and load limits 


— The proposed algorithm provides the best result in terms of profit within 20 seconds. Here the computa- 
tional time is more because the algorithm has been tested for 1000 generations 


5.2. Case 02 
BFO algorithm is tested on the IEEE 30 bus system. In this case, the values of Qo is 300 MW and K 
is 5. During the execution of the algorithm, the code is run for 20 times. Bidding coefficients, output powers 


and profits of the case study are presented in Table 4. It is observed from Table 4 that. 
— All output powers are within generator and load limits 


— The proposed algorithm provides the best result in terms of profit within 20 seconds. Here the computa- 


tional time is more because the algorithm has been tested for 100 generations. 
Error at each iteration for 300 MW and K value of 5 is shown in Figure 4. 


Table 4. Simulations results of IEEE 30 bus system for case 02 
Supplier a b Powers Profit 


1 2.00 0.035525267 128.5317854 524.940744 
2 1.75 0.132519844 36.3426773 151.9170823 
3 1.00 0.208752673 26.66373495 103.979034 
4 3.25 0.081484179 40.69656181 121.1421314 
5 3.00 0.175007819 20.37695205 62.28627279 
6 3.00 0.24496562 14.55765887 46.61630903 
MCP 6.566125934 
Profit 1010.881574 
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Figure 4. Error at each iteration for 300 MW and K value of 5 
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5.3. Comparison and discussions 

The results of bidding coefficients of all methods are tabulated in Table 5. Profits of the proposed 
algorithm are compared with the existing methods and provided in Table 5. Comparison of MCP of all methods 
is given in Figure 5. 

It has been found from these tables that the BFOA yields the best solution. Also, it is observed that the 
computational time is less. The same algorithm can be tested on large scale system in the real-time operation 
of electricity markets. It is identified that the variation of MCP and overall profit is due to 

— The MCP is depending on bidding coefficients. For any small variation in bidding coefficients, it has 
more impact on the MCP 
— The profit also is changed due to variation in MCP 


5.4. Future scope 

In this article, the proposed algorithm has applied for single side optimal bidding strategy problem to 
maximize the profit of the players. The algorithm can be applied for SSOBS by considering security constraints 
on transmission lines. Also, the proposed algorithm may be applied for the same problem with consideration 
of reserve constraints. 


Table 5. Profits of suppliers 
Generator GA[13] PSO[13] FAGSA [13] Proposed 


1 741.45 772.41 1034.9 1208.351 
2 321.32 340.10 376.38 430.013 

3 119.33 125.06 157.22 279.3588 
4 261.01 280.36 498.47 521.2492 
5 125.56 136.32 275.38 324.0376 
6 125.56 136.22 275.38 275.5284 
Total Profit 1694.23 1790.57 2617.73 3038.538 

12 
Marginal cost price 10.15219458 


10 


6.69 


GA PSO FAGSA Proposed 


Figure 5. Comparision of MCP of all methods 


6. CONCLUSIONS 

Single side optimal bidding strategy (SSOBS) has been successfully solved using the BFO algorithm. 
Simulation results in terms of bidding co-efficients, output powers and profits for various loads are provided. 
Also, the results of the proposed algorithm are compared with the available methods. It is observed that the 
proposed BFO algorithm gives better results. It can noticed in the simulation results that The MCP is depending 
on bidding coefficients. For any small variation in bidding coefficients, it has more impact on the MCP. Also, 
the profit is varied due to the variation in the MCP. The proposed algorithm can be used for various other 
problems in operation of power system due to its inherent capabilities in finding the global solution. Also 
the same algorithm can be implemented for real time electricity markets including smart grid with renewable 
energy sources. 
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